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I. INTRODUCTION 


Objectives 

This project was set to examine the evolution of damage density and 
active zone of a crack layer (CL) under different fatigue loading histories in 
a model material. The effect of damage on the near crack tip stress field was 
studied by employing a newly developed technique of semi-empirical stress 
analysis. 

Achievements 

- On the basis of the observed damage growth, we tested the basic hypothesis 
of the crack layer theory, i.e., a self-similarity of damage distribution in 
process of crack layer propagation. 

- A new semi-empirical crack tip stress analysis, based on experimental meas- 

urements of the discontinuity density and the double layer potential 
technique, has been developed as a solution of crack-damage interaction 
problem. Evaluation of the stress intensity factor illustrates the method- 
ology. This illustrative example suggests an alternative to the 

Dugdale-Barenblatt model. 

- The first step in construction of the CL constitutive relationship has been 
made. 

Analysis of the experimental results showed that Arrhenius type constitutive 
relationship describes very well the expansion of the active zone. 

In three consecutive parts we describe our findings. This work is being 
prepared for publication: 


1 


(a) J. Botsis and B. Kunin, 'On Self-Similarity of Crack Layer,' Intern. J. 
Fracture, to appear. 

(b) A. Chudnovsky and M. Ben Onenzdon, 'Semi-Empirical Crack Tip Analysis,' 
Intern, J. Fracture, to appear. 

(c) J. Botsis and A. Chudnovsky, 'On the Expansion of Active Zone of a Crack 
Layer,' to be published. 
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I. DAMAGE ANALYSIS OF A RECTILINEAR CRACK LATER 


Early results of investigations in this program showed that fracture pro- 
pagates as a crack surrounded by a layer of damage, namely, a CL. A well 
developed CL in our model material is shown in Fig. la. 

The self-similarity hypothesis states that the value of the damage den- 
sity p at a point x of the active zone at a time t coincides with that 
at the corresponding point in the initial (t = 0) configuration of the active 
zone, the correspondence being given by a time-dependent affine transformation 
of the space variables 


p t (x) = p 0 (A t * (x - & t )) , (1.1) 

where p^(x) is the value of the damage parameter at the point x at the 

time t, P Q = Pt| . At is a time-dependent 2x2 matrix for a plane problem 

t “0 

(At|^_ being the identity matrix), and the position of the crack tip 

at the time t (JLtlt=0 = °) • Note that » as an y matrix, A t can be uniquely 
decomposed as a product of three time-dependent matrices: a scalar matrix 

(expansion of the active zone), a symmetric positive definite matrix whose 
determinant equals one (distortion of the active zone), and an orthogonal 
matrix (rotation of the active zone). The above formulation of 'self-similar- 
ity' is contained in [1.1] in its infinitesimal form. 

The requirement (1) has two major implications, i) It reduces the evo- 
lution of the function p(x) to the evolution of six scalar parameters: four 

components of the matrix A^. and two components of the vector ii) It 

allows to express the energy release rates associated with the translation 
rotation, and expansion of the active zone (the 'driving forces' in the termi- 
nology of [1.1]) in the form of the well-known J, L an M integrals. 
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respectively [1.4] . The corresponding energy release rate for active zone 
distortion has a similar integral representation [1.1,2], 

A complete examination of the self-similarity hypothesis would require 
comparison of damage distribution within the active zone at different 
instances of its evolution. However, such complete data is not available at 
present. Instead, available is damage distribution within a well developed 
CL. 

Notice that, at any instance t, the trailing edge is the only portion 
of the active zone along which the damage parameter will remain unaltered as 
tne CL continues to evolve. 

For this reason we examine the implications of the self-similarity 
hypothesis for the evolution of damage distribution along the trailing edge. 
The latter is approximated by a straight segment perpendicular to the crack 
path. 

Due to the symmetry of the loading-specimen geometry and the homogeneity 
of the material, (a) the crack follows the straight path along along Xj-axis, 

(b) there is no rotation of the active zone, hence A t is a symmetric matrix, 

(c) the principal axes of A^. coincide with the coordinate axes xj, X2 

(Fig. 1.1). Thus at a time t the position of the crack tip, the gen- 

eric point x of the trailing edge and the matrix A t have the form 



Taking (1.2) into account, equation (1.1) for the points of the trailing edge 
reduces to 

P t (JHt),x 2 ) = p o (0,a 2 (t)x 2 ) . (1.3) 
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(a) General view of a CL in polystyrene, (b) Higher magni 
fication of the boxed area B after thinning. 
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Fig. 1.2 Histograms of craze distribution for eight cross sections 
within the inert zone of a CL. 
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Equation (1.3) implies that, for any two positions of the trailing edge, 
the corresponding along-the-trailing edge damage distributions are related by 
scaling (in x 2 -direction) . 

This restricted form of the self-similarity hypothesis was tested using 
data on crazing distribution in a CL grown under fatigue in polystyrene load- 
ing [1.6] . Quantitative damage analysis was carried out on a thinned specimen 
employing optical microscopy and principles of quantitative stereology [1.5]. 
A micrograph of a typical portion of the CL is shown in Fig. T.lb. Details on 
the experimental procedures can be found in [1.6]. 

Figure 1.2 represents the histograms of craze distribution in vertical 
cross sections at eight locations along the crack path. As a result of the 
present symmetry, the odd central moments of damage distributions in the ver- 
tical cross sections are negligibly small. As simplest indicators of whether 
the shape of the craze distribution stays the same up to a unifo r m dilation 
from section to section, we take the following two: a) the ratio of the total 
amount of crazes for a cross section to six times the standard deviation of 
the damage distribution for the same cross section (this is, for practical 
purposes, the average damage density in the cross section); b) the ratio of 
the variance of the cross-sectional damage distribution to the square root of 
its fourth central moment. Both ratios are independent of scaling and thus 
should remain constant, if 'self-similarity' holds. 

Figure 1.3 (1.4) represents the values of the first (second) of the above 
ratios for the eight aforementioned sections. From the data in Figs. 1.3 and 
1.4, we conclude that both ratios stay constant to within experimental error. 
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Fig. 1.3 For each of the eight cross-sectional craze distributions 
from Fig. 1.2, the 'average craze density' is shown; here 
E is the total number of crazes in a cross section, 
is the standard deviation of the craze distribution in the 
same cross section. 



Crack Length, mm 

Fig. 1.4 For each of the distributions (Fig. 1.2), the ratio of its 
variance y2 t0 the square root of its fourth central 
moment is shown. 
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II. SEMI-EXPERIMENTAL STRESS ANALYSIS 


1. Introduction 

Determination of near and far fields for a crack surrounded by an array 
of microcracks is a key problem for modelling of crack propagation and stabil- 
ity. We approach this problem in two steps: at first a multiple crack 
interaction is characterized experimentally in terms of Crack Opening 
Displacements (CODs). Then, as a second step, we reconstruct the displace- 
ment, strain and stress fields by analytical means employing the results of 
the first step. 

The problem of multiple crack interaction has been recently addressed by 
various authors (a brief review can be found in [XI. 1—3] ) • This problem can 
be expressed formally as a system of singular integral equations with respect 
to unknown CODs. These equations represent the boundary conditions on the 
surfaces of the cracks. Constructive solutions of the problem have been for- 
mulated for just a few types of crack configurations. In general, for a 
random configuration of a large number of microcracks the solution implies 
extremely tedious and time-consuming numerical procedure [II. 4-6]. 

In this section, we employ experimentally observed COD'S as the solution 
of the multiple crack interaction problem. Then the near and far fields will 
be reconstructed by means of the double layer potential technique. The foun- 
dation of the analysis is outlined in the second section and the evaluation of 
an effective Stress Intensity Factor (SIF) reflecting crack-microcrack inter- 
action is presented in the third section as an illustrative example. 
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2. Formulation of a Seml-Bapirloal Analysis 


Let us consider a microcrack of length 2£ at a distance 5 from the 
main crack tip (Fig. II. 1). In this work, we assume plane stress; the case of 



Fig. II. 1 A single microcrack near a main crack, illustrating the 
coordinate system employed. 

plane strain can be treated similarly. The microcrack is characterized as a 
discontinuity bm , i.e., COD (double layer potential). As a result, it 
induces an additional opening displacement on the main crack. Adopting the 
symbols of the previous work of Chudnovsky, et al. [II. 6-8] , we express the 
displacement field u( x ) generated by a microcrack as 


ijmic ( ^ = h £($.x) d? 


(II. 1) 


where oi Is the microcrack line (surface in 3-D) and $(£,x) is the second 
Green's tensor which is defined as the displacement response at the point x 
due to a unit discontinuity at the point of discontinuity % (Fig. II. 1). 
For example, in plane stress, the second Green's tensor for an infinite plane 
is given by (see, e.g., [1.8]) 


<M$.x) = - 


(1+v) 

4nR 2 


1-v 

(n«R - Rn* + n*RI) + 2 


n^R 


RR 


(II. 2) 
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where is the unit normal vector to the line (surface) across which the 
discontinuity takes place, v is Poisson's ratio, I is the unit second rank 
tensor and R is the position vector, i.e., J = £-x (Fig. II. 1). 

The displacement u^ at an arbitrary point x caused by a Microcrack 
Array (MA) is the sum of the displacements generated by each microcrack, i.e.. 


u«(x) 

ivA 


- £ 

k=l 


(k) 

1 imic^) 


■ if* 

k=l /(k) 


(k) (?)<K5,x) d$ 


(II. 3) 


where N is the total number of microcracks in the array. If the size of the 
microcracks is small in comparison with the domain of the array and the number 
N of microcracks is sufficiently large, the summation in Eq. (II. 3) can be 
substituted by integration of an appropriately introduced microcrack density 
over the array domain. Let us introduce a mesh within the domain (Fig. II. 2a) 
and consider a typical square (with an area A Q p = A5 la A5 2 p) which contains 
N Q p microcracks (Fig. II. 2b). Then the displacement vector due to the MA can 
be written as 

») b) 



T 


OCX 

N 

<3 


i 


Fig. II. 2a Schematic representation of the Microcrack Array (MA) surrounding 

the main crack and the subdivision of the domain into a rectili- 
near mesh. 


Fig. II. 2b A typical square for the determination of microcrack opening 
density. 
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N. N* N 


1 ”2 4 ap f 

u A (x) « £ £ £ / b (k> <t)*(?.*> d? 

0=1 P=1 k=l J 


(II. 4) 


- (k) CA ap 


where and N 2 are the numbers of rows and columns in the mesh, respec 

tively. For infinitesimal squares, l{-xl is much larger than the maximum of 

o o 

|A$I, then f(£,x) can be approximated by $U lo , S2p» <n> c p. & * in which 
5j a and (Jjp are the coordinates of the center of the square and <n> a p is 
an average orientation of the microcracks within the square (Fig. II. 2b) . 
Then using the mean value theorem, Eq. (II. 4) can be rewritten as 


N 1 N 2 


^ S^lo'^pWlo'^p'^oP'S^lo^p 
a=l P=1 


(II. 5) 


where the vector £ represents the microorack opening density 


<■£> ■ <A ‘i. 4 V 1 ^ / S ,k>< y 


u (k) C A ( 


op 


(II. 6) 


This can be measured directly as the ratio of the area of the opened cracks 
within a square (shadowed area in Fig. II. 2b) and the area of the square. The 
sum in Eq. (II. 5) becomes an integral over the entire domain of the MA after 
an obvious limiting procedure (|A$| ->-0, N -►•) : 


u * ( x) = / c(£)<H£.<n>,x) d£ 

Mil Af fV <Sl 


(II. 7) 


where V A represents the volume of the MA domain. The corresponding stress 
field <* A (x) induced by the mierocracks array interacting with the main crack 
is presented as 


Za<! } - lx 


/ c(£) <!>(£, <n>,x) d£ 

J 


(II. 8) 
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where J, x i S the stress operator which transforms the displacement field 
into stress 

= T ijk u k I T ijk = ^ (8 ik a j + 8 jk a i ) + X8 ij a k (II. 9) 

Here p and X are Lame's constants, 8jj is the Kronecker's delta (i,j = 
1,2) and k = 1,2,3. The subscript x in (II. 8) indicates that the differen- 
tiation in (II. 9) is to be performed with respect to x. 

3. An Illustrative Rramnl ft; Evaluation of the SIF 

In this section we evaluate the SIF at the crack tip caused by combina- 
tion of the externally applied stress and the MA induced stress <r A . This 

particular example is suggested and closely related to the experimental obser- 
vation of an array of crazes preceding and surrounding fatigue crack in 
polystyrene (Fig. II. 3) [II. 9] . 



Fig. II. 3 An optical micrograph displaying the Crazed Zone (CZ) surrounding 
the crack tip in PS. 

Choosing a cartesian coordinate system with the origin at the crack tip 
(Fig. II. 1) , one can express the SIF as follows 

. r° (ii. io) 

K eff = / G(xj)p(xj) dxj 

/nL _ L 
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Here L is the length of the main crack; GfXj) is the Green's function for 
the SIF due to a dipole unit force applied at x^. In case of a crack in a 
single edge notch (SEN) specimen under mode I loading, G(x 1 ) is given in 
[II. 101. The term p(Xj) is the traction at the line of the main crack 
caused by the external load in combination with the MA interacting with the 
crack, i.e., 

p(Xj) = [a^ + (II. 11) 


where n is the unit normal vector to the main crack face. Thus the total 
(effective) SIF is naturally decomposed into a sum of the conventional SIF K 0 
due to the externally applied load and the SIF K A due to the MA, i.e.. 


- K 0 ♦ K 4 

where 

i r° 

k a = J G<*1>£a ( £i>2x dx l 


( 11 . 12 ) 


(11.13) 


is the stress due to the entire MA expressed by Eq. (II. 8). The SIF K A 
is then given by 



G(xj) 



° ( i )F( i' < 2? >,x l ) d £ dx l 


(11.14) 


c(£) is the microcrack opening density defined by Eq. (II. 6) and F(£,n(£) ,Xj) 
is T x [<K ? ) # <n t >.x]n x at X2 = 0. 

In the case considered, i.e., the crazes are parallel to the main crack, 
<J3^> is the unit vector perpendicular to the main crack for all points 
^1»52>* Then the function F({,<n^>,xi) in E< 1* (11.14) can be written as 
(see Appendix) 
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(11.15) 


F ( £ i » ?2 * ) 


E 

r1 

X 

\ 

rt 

1 

4 7 2 4 

) 4 + 6(£j - Xl ) 2 £ 2 - 3£ 2 

4n 


- Xl > 2 + 

3 J 


where E and v are Young's modulus and Poisson's ratio, respectively. In 
addition, we assume that the MA size at the vicinity of the main crack tip is 
small with respect to the main crack length as well as with the distances from 
the specimen edges. Hence the second Green's tensor for an infinite plane 
(Eq. (II. 2)) can be employed to approximate the second Green's tensor for the 
considered boundary value problem. 

It is useful to introduce the Green's function GgiF for the SIF due to a 
unit discontinuity at the point (£ 1# £ 2 ) normal to a line parallel to the main 
crack. Using Eqs. (11.14) and (11.15) and dimensionless coordinates x 2 , £j 
and $2 representing xj, £1 and £ 2 normalized by the crack length L, 
respectively, we obtain 


9 t 

G SIF*$i»?2* 



G(Xj)F(£ 2 ,£ 2 ,x 2 ) dJCi 


(11.16) 


Depending on the location of the discontinuity (£j,£ 2 ) G SIF can be positive 
or negative. A positive Ggjp reflects an amplifying effect of the discontinu- 
ity, while a negative Ggjp indicates shielding. A plot of the contours of 
equal values of the Ggjp near the crack tip is shown in Fig. II. 4. The con- 
tours are calculated using a numerical integration of Eq. (16) for L/B = 0.5. 
We note that this plot is very similar to the one obtained by Shiue and Lee 
C II. 1 ] using the energy method for the case of a normal strain dislocation 
dipole. Pose [II. 2] and Rubinstein [II. 11] also reported similar results. 
The effect of a microcrack is symmetric with respect to the crack line. The 
border line between the amplification and shielding is a curve which 
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approaches the origin at an angle of about 69°. Apparently, the amplification 



Fig. I I. 4 Contours of equal level of Green's function GgiF for the bIF 

due to a unit discontinuity, normalized by 10^ E/irL/b 0 u). and 
^2 are the coordinates of the discontinuity, normalized by the 
length of the crack L. Here L/B = 0.5. 

zone for mode I loading is located ahead of the crack tip. 

In the considered case of an array of horizontal crazes [II. 9] , the only 
detectable non-zero component of the vector of craze opening density c(£) is 
the vertical component Cjt^) , which is a product of craze opening b 2 (£) 
and the craze density p(£) . The latter is the total length of the craze mid- 
dleline per unit area. This parameter has been studied by Botsis, et al., by 
means of optical microscopy [II. 9]. An example of equal craze density con- 
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tours taken from [II. 9] is shown in Fig. II. 5. Using these data we employ an 
exponential -type approximation of the craze distribution 



Fig. II. 5 


2 3 

Contours of equal craze density in mm /mm [II. 9] for crack length 
L = 1 .9 mm. 


where & a and w are the length and the width of the crazing zone (CZ) ahead 
of the crack tip, respectively. p 0 is the maximum craze density at the crack 
tip estimated as 1400 mm^/mm^ and b Q is the average of the craze openings 
b 2U) • The coefficients a^ are evaluated from the craze density measure- 
ment [II. 9] and found to be 



The overall effect of the CZ can be readily evaluated employing G SIF of 
Eq. (11.16) and the measured craze opening density c(jp approximated by Eq. 
(11.17) 
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(11.18) 


k a = J g sif£ ( 1 ) 

V A 

In the example considered, the following parameters are taken from 
[11.91: crack length = 1.9 mm in a SEN specimen of 20 mm width. Young's modu- 

lus E = 2.2 GPa, Poisson's ratio v = 0.3, and the applied stress a a is equal 
to 16 MPa. In this case K 0 = 1,378 kN/m^ 2 and K A is found to be 
-1128x10® b Q kN/m 5 ^ 2 , in which b Q stands for an average craze opening in m 
upon load application. Since the latter was not measured, we use it as an 
adjustable parameter. Apparently, the total SIF K ef ^ vanishes when b Q = 
1.22x10“® m. The value of 1 to 2x10”® m for craze openings under load seems 
quite realistic, thus suggesting an alternative to the well-known 
Dugdale-Barenblatt proposition [11.12,13], i.e., K 6 *’*’ = 0 due to craze zone 

surrounding the crack. 

5. Conclusion 

The proposed crack tip analysis is essentially based on experimental 
measurements of the discontinuities and/or a discontinuity density. Then 
parameters such as SIF, energy momentum tensor and energy release rates can be 
evaluated by straightforward calculations based on the double layer potential 
technique. 

The semi-empirical method can be also employed for stress analysis of 
crack damage interaction in non-polymeric (non-transparent) materials. With 
this in mind, various polymers and their blends may be selected to serve two 
purposes: a) mimicking brittle as well as ductile response to stress concen- 

tration at the crack tip and and b) to allow an observation and measurements 
of the discontinuity. 
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APPENDIX 


The displacement generated by a microcrack is 

^(x) * y*bj< 


d5 


(Al) 


In case of an horizontal microcrack b 1 (5) =0 and therefore 


u l ( £> = y* b 1 


u 2 (x) 


U>$ 12 (?.x) d 5 


b 2 ( * )$ 22 ( £' x) d 5 


a) 


For plane stress, the stress-displacement relation is given by 


*22< x) = (^2) [u 2,2 ( £> + vu l,l (x)1 


( A2) 


(A3) 


Therefore, the stress operator applied to the second Green's tensor <J>(£,x) 
gives 


T y [$>($,x)] = — 

‘ ~ ~ ~ (l-v 2 ) 




( A4) 


From Eq. (2) in the text, we write the second Green's tensor as 


* 12 ( 1 *£> = ' 


(1+v) 


4nR i 


R 1 R 2 


l-v 

R, + 2 

1+V 1 d 2 


$ 22 ( l' x) = 


(1+V) 

4nR 2 


1-V n 2 

R 2 + 2 — 

R 2 


1 + V 


where Rj = - Xl . R 2 = ? 2 - x 2 and R 


2 = 


( A5) 


2 2 

R 1 + R 2 The derivatives with 


respect to x^ and to x 2 gives 
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4 > 




4, 22,2 ( |'i ) “ 


( l+v) 
4nR 6 


(1+v) 

4nR 6 


[S (•:-«• ssH] 

[&(*••«- 5M] 


Therefore, Eq. (A4) becomes, for x 2 = 0, as 


E 4 2 2 4 

= 7 ( Ri ^ ^ 6R1R9 * 3R9) 

4nR 6 122 


( A6) 
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III. EXPANSION OF THE ACTIVE ZONE 


The requirement of self-similarity hypothesis (Part I) has two implica- 
tions: (i) the evolution of the active zone can be approximated by four 

elementary moments, namely, translation and rotation as a rigid body, expan- 
sion and distortion. (ii) It allows to express the active parts of the 
corresponding driving forces in the form of the J, L, M and integrals, 

respectively [III.l]. 

Crack propagation which coincides with the translation of the active zone 
is a relatively fast process. The CL theory derives the kinetic equation for 
active zone translation by employing the principle of minimum entropy produc- 
tion [111.1,2], On the other hand, active zone deformation and rotation 
result from slow processes of damage growth. 

On the basis of experimental observations, various propositions can be 
made with respect to the constitutive equations for expansion, distortion and 
rotation of the active zone. Herein we concentrate on the expansion of the 
active zone only. 

In this respect we have examined two candidates for constitutive rela- 
tions; Onsanger's type linear relation 

e = LX exp (III.l) 

and Arrhenius type exponential relation 

e = Cj exp [-C 2 X exp ] (III. 2) 

between the rate of expansion e, and the corresponding force X exp . The 
coefficients L[ sec •Joules -1 ] (Eq. (III.l)), Cjfsec] and C 2 t Joules -1 ] (Eq. 
(III. 2)) are obtained from linear regression analysis. 
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Under the assumption of homogeneous deformation within the active zone. 


the rate of expansion e with the crack tip as the origin is 


where w and Jla are the width and length of the active zone, respectively 
(Fig. I.l). The force X exp is expressed as the difference between active 

and resistive parts, i.e., 

X exp = M " 1*0 (III ‘ 3) 

Here M is the energy release rate due to expansion of the active zone, y 
is the specific enthalpy of damage and R 0 is the resistance moment of the 
active zone associated with the expansion. Both M and R Q are approximated 
as 

M = \ oJ,Jta (III. 4) 

2 1 

R 0 = <p 0 >A (III. 5) 

where Jj is the energy release rate due to translation of the active zone, 
£a is the length of the active zone and a is a dimensionless coefficient 
with 0 < ail. <p 0 > is the average damage density within the active zone and 
A its area. 

Analysis of CL stability has shown that the critical energy release rate 
Aj c , can be expressed as the product of the specific enthalpy of damage y * , 
and the resistance moment at critical propagation, namely, A lc = y*Ric* This 
process is a manifestation of mainly new damage nucleation. Thus y* is 
directed related to the corresponding specific energy (i.e., energy of nuclea- 
tion of new damage). On the other hand, y which appears in Eq. (III. 3) 
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corresponds mainly to the energy of relative slow process which follows the 
nucleation. Furthermore, it is well recognized that the difference between 
the energy of new phase nucleation and the energy of phase growth is essential 

[111. 3] . Accordingly, it is expected that y* is significantly greater than 

Y. 

For our model material, y* has been found to be of the order of 30 J/m^ 

[111. 4] . An experimental method to evaluate y is being developed. To com- 
pare the two candidates (Eqs. (III.l) and (III. 2)) for constitutive equations 
for active zone expansion, we take y to be 10% of y*. 

Fracture propagation kinetics are observed under two different fatigue 
load histories. Whereas the frequency and load ratio were the same, the level 
of mean stress is 16.0 MPa and 10.7 MPa, respectively. 

Analysis of the experimental results according to Onsanger relationship 
(Eq. (III.l)) resulted in two different values of the kinetic coefficient L, 
namely, L = 0.98xl0 -2 sec • joules -1 for the first experiment and L = 
0.68x10 - 2 sec •joules -1 for the second experiment. In addition, in both cases 
the coefficient of correlation in the linear regression analysis, in the e, 
*exp Plane, was low, regardless of the value of the parameter a. 

This suggests that simple Onsanger relationship is not adequate. 

The solid lines in Fig. III.l represent the right-hand side of (III. 2). 
The data points are measurements of the rate of expansion of the active zone. 
The correlation coefficient for both set of data was of the order of 0.9 and 
the pre-exponent parameter Cj was 1.12xl0 -3 sec and 0.34xl0 -3 sec, respec- 
tively. On the other hand, the coefficient C 2 was found to be the same for 
both experiments regardless of the parameter a. If an Arrhenius type consti- 
tutive relationship was to describe the expansion rate, C 2 should be the 
same in this case since both experiments were performed under different load 
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levels and the same temperature. This Is very strong evidence that Eq. 

(III. 2) could be adopted to describe the expansion rate of the active zone. 
In order, however, to unquestionably assess the applicability of an Arrhenius 
type kinetic equation for active zone expansion, experiments with different 
loading rates and temperatures should be available. 



Fig. m.i Rate of expansion as a function of the expansional thermodynamic 
force for two levels of mean stress. Note that the lines are 
parallel. 
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